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Roton-Rotation Coupling of Acetylene in 4 He 

R. E. Zillich and K. B. Whaley 

Department of Chemistry, University of California, Berkeley, CA 94720 

Rotational absorption spectra of acetylene in superfluid 4 He are calculated using a combined 
correlated basis function-diffusion Monte Carlo method (CBF/DMC) as well as a numerically exact 
path integral correlation function approach (PICF). Both approaches predict anomalously large 
distortion constants in addition to reduced rotational constants, and the higher rotational excitations 
are seen to be accompanied by an additional peak and absorption band. The CBF analysis shows 
that these phenomena are due to strong coupling of the higher rotational states with the roton and 
maxon excitations of 4 He, and the assocated divergence of the helium density of states in this region. 

PACS numbers: 05.30.Jp, 33.20.Bx, 33.20Sn, 36.40.-c, 67.40.Yv 



Helium nanodroplets have become a useful matrix for 
isolation spectroscopy Unlike in classical liquids, the 
spectra of embedded molecules may be rotationally re- 
solved in the superfluid environment of 4 He, which has 
led to numerous investigations of rotational constants of 
molecules in 4 He droplets (see table 1 in Ref. 0). For 
all but the lightest rotors, a significant reduction of the 
rotational constants with respect to their gas phase val- 
ues has been found, with identical symmetry, which has 
been explained by the adiabatic following of a molecular- 
interaction induced non-superfluid fraction of the solva- 
tion shell. In contrast, the lightest molecules exhibit 
relatively small reduction of their rotational constants 
in 4 He droplets Calculations have shown that the 
assumption of adiabatic following breaks down for light 
molecules 4] and that instead of coupling to localized 
helium solvation density, there is appreciable coupling to 
delocalized collective helium modes [fg. 

Acetylene (HCCH) is one such light molecule, whose 
ro-vibrational spectrum has been obtained by infra-red 
spectroscopy in 4 He droplets [|| . The effective rotational 
constant B e g was found to be reduced to 88% of its gas 
phase value (1.1766 cm -1 ) and a very large effective dis- 
tortion constant D e g was measured, showing a relative 
enhancement D e f{/B an order of magnitude larger than 
that seen for heavier molecules such as OCS, SFg and 
(HCN)2. It was speculated that the excessive broaden- 
ing of the IR spectral lines might indicate coupling to 
excitations of the helium, but in the absence of theoreti- 
cal understanding of such a coupling no conclusions could 
be drawn. The calculations of the rotational absorption 
spectra of HCCH presented here show how coupling be- 
tween the molecular rotation and the roton excitations of 
4 He results in an anomalously large distortion constant as 
well as a reduced rotational constant, with predicted val- 
ues in good agreement with experimental measurements. 
They also reveal that the higher rotational excitations 
display a secondary peak as well as a weak maxon-roton 
band between them that both result from this coupling. 
We discuss how these additional spectral features might 
be observed, in order to stimulate experiments to mea- 
sure this unique signature of rotation-roton coupling. 



We present two independent calculations. The first 
calculation is a combination of correlated basis func- 
tion (CBF) theory and diffusion Monte Carlo (DMC) 
methodology that derives zero temperature excitation en- 
ergies from linear response theory for correlated wave 
functions and then evaluates these with numer- 
ically exact DMC calculations for the required in- 
put ground state properties |5j. The rotational ex- 
citations are described in this CBF/DMC approach 
within an analytic approximation that reveals the mi- 
croscopic origins of the rotational excitations and as- 
sociated spectral absorption features. For the second 
calculation we use path integral Monte Carlo meth- 
ods H, (at temperature T = 0.625 K) to evaluate 
the imaginary-time orientational correlation functions 
Fj(t) = (l/Z)J2 m Tr{Y/ m m))Yj m (n(0))e- H / kT } (fi 
is the molecular orientation), sampled from the full 
molecule-helium Hamiltonian H that includes molecular 
rotation and all translational degrees of freedom as well 
as the permutation exchange symmetry of helium |lfj| . 
The inverse Laplace transform (made with the maxi- 
mum entropy method 01) °f ^j(t) is proportional to 
the rotational absorption spectrum for dipole (J = 1), 
quadrupole (J = 2), and higher order transitions. The 
PICF calculation provides a quantum simulation free 
from systematic approximation which is important for 
verification of the novel spectroscopic features predicted 
by the CBF/DMC calculations. 

In Ref. 5], we formulated the CBF/DMC approach 
for molecular excitations in 4 He and showed how CBF 
theory is able to describe the effect of coupling of collec- 
tive excitations of the 4 He environment to the rotational 
dynamics of a molecule. This coupling was seen to renor- 
malize the rotational energy Tiujj corresponding to total 
angular quantum number J for a linear molecule with 
bare mass M and gas phase rotational constant B by an 
energy-dependent self energy Ej(o;j) [l2| . 



Tiuj ; 



BJ(J + 1) + 



(1) 



which is a functional of several ground state quanti- 
ties The Legendre coefficients ge>(r), £' > 0, of 
the pair distribution function between HCCH and 4 He 



2 



are obtained from a DMC simulation of HCCH and 256 
4 He atoms (periodic boundary conditions), employing a 
trial wave function of the form used in Ref. |5J (parame- 
ters c = 7.724 A and b = 2.670 A), interaction potentials 
from Refs. 0, 0] , and implementing descendent weight- 
ing 0. Odd components of gi>{r) vanish for HCCH 
and the zero component gt=o{r) does not contribute to 
Sj(w), so that the £' = 2 coefficient is the dominant con- 
tribution to Tij(uj). From angular momentum conserva- 
tion [j| , one can then show that for HCCH the dominant 
effect of the superfluid 4 He environment on the molecu- 
lar rotational excitations J is the coupling of molecular 
rotational states that differ by | J — £\ = 2. 

In general, the self energy is complex, leading to finite 
life-times and decay of rotational excitations into other 
modes. When the momentum p is such that the input 
energy tko equals the sum of three excitation energies 
B£(£ + 1) (molecular rotation), e(p) (phonon-roton ex- 
citation), and h 2 p 2 /2M (kinetic energy of the recoiling 
HCCH), these states can be excited by decay of the ro- 
tational excitation hu> = fitoj. 

Eq. J3> needs to be solved self-consistently for ujj since 
the self energy Ej(wj) is itself a function of ujj. When 
the imaginary part of Sj is small, the real part of the 
solutions are seen to correspond to the positions of well- 
defined peaks in the dipole (J = 1), quadrupole (J = 
2),. . . absorption spectra of the molecule from its ground 
state. These spectra, normalized to unity, are given by 



Sj(uj) = -$m[BJ(J+l) + Ej(w)-fiw]- 
9;m£j(o;)/7r 



( 7 j(w) - M 2 + (StoEjM) 5 



(2) 



where we have defined r y,j(uj) = BJ(J + 1) + 5Je£j(w). 
The imaginary part of the self energy, SmSj, gives 
the homogeneous width of the Lorentzian line associated 
with decay of the excitation of life-time fi/3mSj(uj) 
into the lower energy modes of the molecule (rotation 
and translation) accompanied by emission of a 4 He exci- 
tation. In the opposite limit when the imaginary part of 
Sj(w) is large, eq. (J2J) shows that the absorption Sj{lo) 
can become large regardless of whether or not uj is a so- 
lution of eq. CJ. In this case we find a broad band which 
cannot be uniquely associated with a single molecular ro- 
tation mode, but instead is due either to a simultaneous 
excitation of both the molecular rotation and 4 He modes, 
or of only 4 He modes (see below). 

Fig. ^ shows the ground state absorption spectra of 
HCCH for J = 1,2,3 calculated by CBF/DMC, eq. (2). 
The rotational excitations obtained from PICF are also 
shown, as broad solid lines with height proportional to 
their spectral weight. Lineshapes and broad spectral fea- 
tures cannot be obtained with the correlation function 
spectral inversion approach, but the positions and rela- 
tive strength of discrete peaks are very reliable 01 ■ The 
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FIG. 1: CBF/DMC absorption spectra Sj(u), J = 1,2,3, 
for HCCH in bulk 4 He (thin solid lines) and PICF rotational 
excitations (broad solid lines, height proportional to the spec- 
tral weight) for HCCH in a cluster of 64 4 He. The CBF/DMC 
spectra have been broadened by 10 mK. The bottom panel 
shows the density of states e'(p) + TL 2 p/M. The vertical dotted 
lines in the J = 2 and J = 3 panels indicate the edges of the 
roton-maxon bands coupling these states to I — and I = 1 
states, respectively. 



CBF/DMC spectra have sharp Lorentzian peaks associ- 
ated with excitations of finite life-time, as well as broad 
features between these peaks. The peaks are in good 
qualitative agreement with the PICF excitations, which 
are made for HCCH in a cluster of 64 4 He atoms (ap- 
proximately two solvation layers). This implies that the 
relevant 4 He collective excitations are already well devel- 
oped at this cluster size and that the rotational coupling 
to these modes is not sensitive to temperature variation 
below T ~ 0.6 K. 

We analyze first the sharp peak structures. The loca- 
tions Tito' of peaks in the CBF/DMC spectrum Sj(u>) are 
found by solving the equation Jj(ui) — hto = (see also 
eq. EJ, which can have one or more solutions. Each solu- 
tion fiu)' can be associated with a rotational excitation of 
finite life-time, which can decay into a £ < J state and an 
excitation of the helium environment. It is evident from 
Fig. that while for J = 1 a single peak is obtained, 
the CBF/DMC absorption spectra Sj(oj) for J = 2 and 
J = 3 have two peaks each, in contrast to the expected 
pattern of excitations for a linear rotor which has only a 
single excitation at energy BJ(J + 1) for a given J. This 
surprising double peak feature is also found for the rota- 
tional excitations calculated with the numerically exact 
PICF method. The latter do not rely on any analytic 
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FIG. 2: The absorption spectrum Sj(u>), for the J = — > 2 
transition of HCCH in 4 He, where we have varied the gas 
phase rotational constant B between 0.5 x B (bottom curve) 
and 1.5 x B (top curve), in steps of 0.1 x B. 
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TABLE I: Rotational excitation energies of HCCH in 4 He 
calculated with CBF (bulk 4 He) and PICF ( 4 He 64 ), for J = 
1, 2, 3, compared with experimental values extracted from IR 
spectra J' denotes the secondary peaks and w the respec- 
tive spectral weights. 



approximations and thus provide additional verification 
of the correctness of the splitting seen in the CBF/DMC 
spectra. We therefore are led to conclude that unlike 
heavier rotors such as OCS, the pure rotational spec- 
trum of HCCH in 4 He for excitations higher than J = 1 
does not conform with the simple pattern predicted for 
an effective linear rotor. 

Inspection of the CBF equations allows us to deter- 
mine that the reason for the occurrence of two peaks for 
absorptions to J = 2 and J = 3 is the non-monotonic be- 
havior of the collective excitation spectrum in 4 He. The 
divergent density of states (see bottom panel of Fig. 
at both the roton minimum (wave number p r ) and the 
maxon maximum (wave number p rn ) of the helium excita- 
tion spectrum 0] results in divergences in the self energy 
Ej(lj) at to = B£(£ + 1) + e(p) + h 2 p 2 /2M for p = p r ,p m 
and I — 0,1,2,.... These divergences can split spectral 
peaks that are close to the roton and maxon into two, 
shifting the "primary" peak below the roton minimum 
and moving the "secondary" peak above the maxon max- 
imum The dipole spectrum, J = 1, has only a single 
peak, because the J = 1 excitation energy is much lower 
than the roton energy and consequently does not couple 
so effectively to the collective modes. Consequently the 
secondary J = 1 absorption peak has negligible weight 
and is not observed. For J = 2, the weight of the sec- 
ondary peak is small but non-negligible, while for J = 3 
more than half of the spectral weight is actually carried 
by the secondary peak, see Table [I] 

In order to further understand the coupling between 
rotations and rotons that is implied by this splitting, we 
have calculated the CBF /DMC spectrum as a function of 
the molecular gas phase rotational constant B. This al- 
lows the molecular rotational energy to be tuned across 
the roton/maxon energy regime, and the coupling be- 
tween molecule and helium to be modulated as a result. 
Fig. shows the quadrupole spectrum Sj(uj), J — 2, as 
B is artificially varied between 0.5 x B and 1.5 x B in 
eq. J2J and in the expression for the self energy. For sim- 
plicity we use the same molecule-helium pair distribution 



function for all B values in this comparison. When the 
gas phase value of B is sufficiently small that the J = 2 
excitation does not couple to the roton excitations it is 
evident that Sa(w) now shows only a single peak, just like 
Si(uj) in Fig. n On increasing B to its actual gas phase 
value, 1.1766cm _1 , the secondary peak is seen to emerge 
and simultaneously a weak roton-maxon band starts to 
grow in between the two peaks. On increasing B beyond 
the HCCH value, the secondary peak continues to grow 
until it eventually accounts for the entire spectral weight 
and the primary peak vanishes. 

Between these two peaks in the J = 2, 3 spectra, there 
are no sharp rotational peaks because excitations here 
would decay instantaneously due to the high density of 
states between the roton and maxon excitations of 4 He. 
The resulting large imaginary part of Sj(w) leads to 
a broad roton-maxon band in the spectrum. (Unfortu- 
nately such broad features cannot be retrieved from an- 
alytic continuation of the PICF to real frequencies.) Al- 
though these roton-maxon bands appear small in Fig. ^ 
they do possess non-negligible weights of 0.06 and 0.08 
in the J = 2 and J = 3 spectra, respectively. The roton- 
maxon band for J = 3 is shifted in energy by 2B with 
respect to the corresponding band for J = 2, because 
for J = 2 the molecular rotation is coupled to a pure 
excitation of the helium environment with the molecule 
in its rotational ground state £ — 0, while for J = 3, 
the rotation is coupled to an I = 1 molecular state. In 
both cases J — £ — 2, which corresponds to the strongest 
contribution to the self energy £j(u;), as noted above. 

We now consider how these CBF and PICF results 
may be observed experimentally. Our ground state spec- 
tra correspond to dipole, quadrupole, etc. rotational ab- 
sorption from the J = ground state, while in Ref. 
rotational transitions J = «-» 1 and J = 1 — > 2 were 
observed (accompanied by excitation of the C-H vibra- 
tional stretches). Table [I] compares the CBF and PICF 
rotational excitations with the experimental values ex- 
tracted in 0. The two excitations J = 1 and J = 2 
that were accessed by the spectroscopic measurements 
are seen to be in excellent agreement with the theoret- 
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CBF 


PICF 


Exp. [6] 


gas phase [6] 




1.075 


1.031 


1.042 


1.172 




0.0175 


0.0179 


0.0195 
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TABLE II: Rotation constant B e g and distortion constant 
D off of HCCH in 4 He. The CBF and PICF results are com- 
pared with the fit made toJ = 0— »1,J = 1— * 2 rovibra- 
tional fine structure transitions in IR spectra 6]. (Units in 
cnr 1 ) 



ical values for J = 1 and for the primary peak of the 
J = 2 spectrum. (The CBF self energy and accuracy of 
results can be further systematically improved 0.) Ta- 
ble [n] compares the spectroscopic constants B e g and D e g 
obtained by fitting these two excitations to the non-rigid 
linear rotor spectrum B cS J(J + 1) — D cti (J(J + l)) 2 . 
The distortion constant relative to B, D e g/B, is very 
large compared to that for heavier molecules 0, [6] . Our 
calculations show that this anomalously large distortion 
constant results from the strong coupling of the J = 2 
excitation with the helium roton and maxon excitations. 
Fig. [21 shows that the spectral shift to lower frequencies 
that is responsible for D e g is greater for lighter molecules, 
reflecting an increased coupling as the molecular exci- 
tation become more resonant with these helium modes 
that derives primarily from their high density of states. 
In Ref. || we demonstrated that this deviation from a 
linear rotor spectrum is much smaller when the molecule 
couples to a simple linear (phonon) excitation spectrum. 

In this CBF calculation the homogeneous linewidth of 
the absorption spectra in Figs.C]and[5J 3mSj(w), results 
from rotational relaxation only. For J = 1, the rotational 
linewidth vanishes due to the symmetry of the HCCH- 
4 He potential (see eq. (4.4) in Ref. @). For J = 2, it is fi- 
nite but still smaller than the experimental IR linewidth, 
which is believed to also have contributions from vibra- 
tional relaxation as well as inhomogeneous broadening 
due to finite size effects 6J. Vibrational relaxation may 
lead to a reduced life-time of the secondary transition 
such that it is too short to be observed. Multi-phonon 
excitations, which become prevalent in 4 He for excita- 
tion energies around 20 K and higher, and are not incor- 
porated in the present CBF analysis, may also broaden 
the theoretical peaks. For HCCH the secondary peak for 
J = 2 is predicted to lie outside the range accessed in 
Ref. and a preliminary search for this at higher IR en- 
ergies has not been successful 0] , possibly for the above 
reasons. Another complication for comparison with the 
IR spectra is that the J = 2 state is accessed by the 
dipole transition J = 1 — > 2, while the CBF calculation 
yields the quadrupole transition J = — > 2, with differ- 
ent transition matrix elements for the HCCH- 4 He cou- 
pling. Measurement of the rotational Raman spectrum 
may allow a direct access to the J — — > 2 transition. 

In summary, calculations of the rotational absorption 



spectra corresponding to dipole, quadrupole, etc. transi- 
tions of HCCH from its rotational ground state in bulk 
4 He, -> J, using the CBF/DMC and PICF approaches 
yield an anomalously large distortion constant D c g that 
is seen to result from a strong coupling of the higher 
molecular rotation states to the roton excitation of he- 
lium. The CBF/DMC spectra for J = 2 and J = 3 
show secondary peaks at higher energies, which are con- 
firmed by the numerically exact PICF results, as well as 
weak roton-maxon bands between the two peaks. We 
expect that these unique signatures of coupling between 
roton modes of helium and rotational modes of embed- 
ded molecules will be found in Raman spectra, as well 
as possibly in microwave and IR spectra for other light 
molecules. 
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